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COMPARISON  OF  RECEIVER  OPERATING  CHARACTERISTICS  OF  FOUR 
DIFFERENT  PROCESSORS  FOR  THE  DETECTION  OF 
RANDOM  SIGNAL  TRANSMISSIONS 

INTRODUCTION 

In  an  effort  to  avoid  interception,  a  random  transmitted  signal  is  sometimes  used  for 
detection  and  communication  purposes.  It  is  presumed  here,  for  the  most  part,  that  the  receiver 
(which  may  or  may  not  be  co-located  with  the  transmitter)  has  precise  knowledge  of  the  particular 
transmitted  random  signal  sequence,  but  not  of  the  medium  attenuation.  The  receiver’s  goal  is  to 
detect  any  target  echo,  if  present,  with  a  maximum  detection  probability,  while  maintaining  a 
small,  fixed  false  alarm  probability. 

Various  processors  are  candidates  for  signal  detection  in  this  scenario,  among  which  are 

•  the  energy  detector, 

•  the  crosscorrelator, 

•  the  normalized  crosscorrelator,  and 

•  the  optimum  detector. 

The  last  processor  and  its  detectability  performance  serve  as  a  standard  of  comparison  for  the  first 
three  practical  processors.  Comparison  of  these  four  processors  is  desired  for  time-bandwidth 
products  ranging  from  small  values  to  large  values  and  for  all  ranges  of  received  signal-to-noise 
ratio. 


Under  the  conditions  to  be  indicated  explicitly  below,  the  characteristic  functions  of  the 
decision  variables  for  all  four  processors  have  been  derived  in  closed  form.  A  numerical 
procedure  utilizing  a  fast  Fourier  transform  (FFT)  then  yields  the  entire  exceedance  distribution 
function,  from  which  the  receiver  operating  characteristics  are  readily  found.  This  exact  approach 
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avoids  the  need  to  resort  to  approximations  such  as  the  deflection  criterion,  with  its  limited 
applicability.  The  misleading  character  of  the  deflection  for  the  four  processors  of  interest  will  be 
discussed  later  in  the  report. 
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OPTIMUM  PROCESSORS  FOR  VARIOUS  SIGNAL  MODELS 


The  transmitted  signal  sequence  is  characterized  by  {sm}  for  1  <  m  <  M,  where  Mis  the 
number  of  signal  samples.  This  sequence,  whether  deterministic  or  random,  is  known  to  the 
receiver.  The  received  or  measured  data  sequence,  under  the  signal-absent  hypothesis  H0  or  the 
signal-present  hypothesis  Hu  respectively,  is 


x 


m 


nm  for  H0 
aSm  +  nm  f°r  #1 


for  1  <  m  <  M , 


where  scaling  a  is  the  medium  attenuation.  This  scaling  a  is  constant  over  the  signal  extent  M; 
however,  the  knowledge  of  a  or  its  statistics  may  vary  with  the  situation.  The  additive  noise  {nm} 
is  stationary  white  Gaussian  with  zero  mean  and  known  variance  a2 .  M can  also  be  interpreted 
as  the  time-bandwidth  product  of  the  signal. 

The  optimum  and  near-optimum  processors  will  be  derived  under  a  variety  of  assumptions 
regarding  the  character  or  knowledge  of  scaling  a.  For  comparison,  the  situation  of  a  determin¬ 
istic  transmitted  signal  sequence  is  addressed  first;  then,  the  alternative  current  situation  of  a 
random  transmitted  signal  is  undertaken.  The  derivations  of  the  performances  of  the  resultant 
processors,  in  terms  of  false  alarm  and  detection  probabilities,  are  considered  in  the  next  section. 


DETERMINISTIC  TRANSMITTED  SIGNAL 

In  this  subsection,  not  only  is  transmitted  sequence  {sm}  known  to  the  receiver,  but  it  has 
no  variability  to  it.  Thus,  its  sample  (transmitted)  energy 

M 


is  a  known  constant.  For  a  given  fixed  value  of  medium  attenuation  a,  the  probability  density 
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functions  (PDFs)  of  observation  X=  (x1( . . . ,  xM)  under  H0  and  H\  are 


p0(x)=  n 

m=  1 


(2«*r 


exp 


,2  > 


2  at) 


and 


(1) 


using  the  knowledge  of  {sm}.  For  this  fixed  value  of  a,  the  likelihood  ratio  (LR)  is 


LRW=^frexp 


r  a  a2  ^ 

TT  C  ~  ry  2  Q  ]  > 
yan  2<y  J 


where 


M 

c-E 

m=\ 


Vm 


is  the  crosscorrelation  of  the  known  signal  (sm)  with  the  received  data  {xm}.  Several  cases  of 
knowledge  of  scaling  a  will  now  be  considered. 


Case  I:  a  =  Known  Constant 

In  this  case,  the  quantity  \a2  Qja1n  is  a  known  constant.  Using  the  monotonicity  of  the 
exp( )  function,  the  LR  test  is 
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where  fixed  threshold  T  absorbs  all  the  known  constants  a,  Q,  and  an. 

Case  II:  a  =  Unknown  Constant 

Now,  ja2  Q/o2n  is  unknown,  but  it  is  constant,  that  is,  it  does  not  vary  from  trial  to  trial. 
Therefore,  the  LR  test  is  again 


where  fixed  threshold  T  absorbs  the  unknown  constant  a  and  the  known  constants  Q  and  a„. 


Case  III:  a  =  Random  Variable 

If  random  variable  (RV)  a  has  known  PDF  pa,  the  PDF  p\(X)  in  equation  (1)  must  be 
averaged  over  variable  a.  The  modified  LR  is  now 


ttw^x  \dupa{u)p^X\a  =  u)  fj  /x  , 
LR(^)  = - ~p{5t) - =  dU  Pa^  CXP' 


r  u  u2  N 


(2) 


If  random  variable  a  is  limited  to  positive  values,  this  LR  is  monotonic  with  quantity  C,  which 
may  be  seen  from 


d 

8C 


oo 

LR(X)  =  \dupa  (u)  exp 


u 


Q) 


Because  the  integrand  is  everywhere  positive,  it  follows  that  the  right-hand  side  is  positive  for  all 
C,  regardless  of  known  constant  <9;  that  is, 


LR(if)  =  l(C;  Q,  a1,,  p.) 
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in  equation  (2)  is  a  monotonically  increasing  function  of  C.  Since  Q,  a2n,  and  the  parameters  of 
PDF  pa  are  constant  from  trial  to  trial,  the  LR  test  is  simply 


C*T, 


where  fixed  threshold  T  includes  Q,  <j2n,  and  the  parameters  of pa. 


Notice  that  the  LR  test  here  is  always  a  comparison  of  crosscorrelator  output  C  with  a 
fixed  threshold,  regardless  of  the  PDF  pa  of  positive  amplitude- scaling  a.  If  a  could  take  on 
negative  values,  this  same  conclusion  need  not  hold  true.  Thus,  for  a  deterministic  signal  {sm}, 
the  optimum  processor  is  always  C  T,  regardless  of  the  knowledge  or  statistics  of  scaling  a  (if 
positive).  This  is  the  standard  comparison  of  the  crosscorrelator  output  with  a  fixed  threshold. 


RANDOM  TRANSMITTED  SIGNAL 

Here,  signal  segment  {sm}  will  be  drawn  from  a  random  sequence,  which  need  not  be 
stationary  or  white.  The  particular  values  (sm)  are  again  known  to  the  receiver,  but  will  vary  from 
trial  to  trial.  Now,  the  signal  sample  (transmitted)  energy 

M 


is  known  to  the  receiver,  but  it  varies  from  trial  to  trial.  For  a  given  fixed  value  of  attenuation  a, 
the  LR  is  again 


LR(X)  =  exp 


(3) 


using  the  knowledge  of  the  particular  samples  {sm}.  Again,  three  cases  of  knowledge  about 
scaling  a  are  of  interest. 
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Case  I:  a  =  Known  Constant 


In  this  case,  the  quantity  \a2  Q/<J2n  is  known  at  the  receiver,  but  it  varies  from  trial  to 
trial.  Therefore,  the  LR  test  can  only  be  simplified  to 

C-iaQ^T, 


where  T  is  a  fixed  threshold.  The  quantity  \aQ  must  be  subtracted  from  C  because  it  varies  with 
each  trial;  only  after  this  subtraction  is  a  comparison  with  a  fixed  threshold  valid  as  the  LR  test. 
The  test  C  <  T  is  no  longer  the  optimum  processor  when  {sm}  is  a  segment  of  a  random  process, 
even  though  the  entire  sequence  {sm}  is  known  to  the  receiver. 

An  alternative  way  to  express  the  LR  test  is 


MM 

S  (xm-asm)  <  T . 

m=  1  m=l 


This  expression  is  a  difference  of  two  energy  estimates:  one  with  the  data  {xm}  as  is,  the  other 
with  the  received  signal  subtracted.  The  leading  term  (by  itself  as  a  processor)  is  called  the  energy 
detector. 


If  the  number  of  samples,  M,  in  the  signal  segment  is  large,  the  quantity  Q,  although 
random,  will  not  vary  much  from  trial  to  trial.  In  this  special  case,  the  term  \aQ  will  be 
substantially  constant,  and  the  LR  test  becomes  approximately  C<T  for  large  M,  where  fixed 
threshold  T  absorbs  the  almost-fixed  term  \aQ . 

Observe  that  the  statistics  of  signal  {sm}  were  never  needed  for  the  above  derivations  of 
the  optimum  processor;  these  results  apply  for  any  segment  of  any  random  sequence  {sm}  known 
to  the  receiver.  However,  the  performance  of  the  LR  test  will  need  specification  of  the  signal 
statistics. 
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Case  II:  a  =  Unknown  Constant 


Here  the  quantity  jaQ  is  unknown  and  varies  from  trial  to  trial.  Therefore,  the 
calculation  of  C  -  jaO  cannot  be  carried  out,  and  the  optimum  processor  cannot  be  realized, 
because  a  is  unknown. 

When  the  signal  samples  {sm}  for  1  <  m  <M  happen  to  have  a  high  level,  then  Q  is  also 
high  level  and  the  subtracted  term  tends  to  compensate  for  that  effect  on  C.  Similarly,  a  low-level 
signal  sample  set  is  subject  to  the  same  type  of  compensation.  That  is,  the  optimum  processor  is 
compensating  for  the  effects  of  a  variable  signal  transmission  level.  This  suggests  the  following 
suboptimum  processor  that  also  compensates,  but  does  not  need  knowledge  of  unknown  gain  a , 
namely, 


C  > 

4q<t’ 


where  T  is  a  fixed  threshold.  The  alternative  equivalent  expression. 


I  with  rm  -  ~j=  for  1  <  m  <  M , 


reveals  that  the  local  reference  {rM}  is  independent  of  the  absolute  level  of  the  particular  signal 
segment,  {sm},  employed. 


Another  approach,  in  this  case  of  an  unknown  parameter,  uses  the  method  of  maximum 
likelihood  (ML).  The  a-dependency  of  PDF  p^X)  is  given  by  equation  (3);  therefore,  the  ML 
estimate  is  a  =  0  for  C  <  0  and 


m=l 


for  C  ^  0 . 
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This  choice  reflects  the  fact  that  a  is  known  to  be  positive.  The  corresponding  generalized 
likelihood  ratio  (GLR)  is 


GLR(X)  = 


exp 


1 


_  £ 

Ucr„2  Q) 

1 


for  C  >  0 

for  C  <  0 


which  is  monotonic  in  C.  The  GLR  leads  to  the  test 

~Jq><T  (wherer>1)> 

in  agreement  with  the  compensation  argument  suggested  above.  This  processor  is  called  the 
normalized  crosscorrelator  or  normalized  matched  filter  for  random  signal  transmissions. 


It  should  be  observed  again  that  knowledge  of  the  signal  statistics  was  not  needed  or  used 
in  arriving  at  suboptimum  processor  Cj yfQ.  Also,  the  earlier  comment  relative  to  a  large  sample- 
size  M  is  relevant  here  as  well;  namely,  Q  is  relatively  constant  then,  and  the  LR  test  is  well 
approximated  by 

C<  T  for  large  M, 
where  T  is  a  fixed  threshold. 


Case  III:  a  =  Random  Variable 


The  LR  is  again  given  by  equation  (2);  that  is, 


LR(X)  =  jdu  pa(u)  exp 


u 


vcr; 


c 
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But  now,  both  C  and  Q  vary  with  each  trial;  it  is  no  longer  sufficient  to  consider  d(LK)/dC. 

In  order  to  proceed  further  in  this  case,  it  is  necessary  to  specify  the  particular  PDF  pa. 
Let  us  consider 


(“-”Q2 

2<7a  _  ’ 

where  mean  ma  and  variance  <x02  of  Gaussian  RV  a  are  presumed  known.  (If  ma  >  3cra,  the  RV  a 
is  almost  always  positive,  which  is  the  situation  of  interest.)  The  LR  then  becomes 


PM  =  (2^a)  12  eXP 


LR(X)  = 


exp 


ml  (o-X  + 

2cr2a  2a2aa2n(a2n  +  a2aQ) 


The  quantity  C  is  more  variable  than  Q  from  trial  to  trial.  Also,  the  exponential 
accentuates  variations  far  more  than  the  square  root  in  the  denominator.  Therefore,  the  LR  test  is 
well  approximated  by 

(°X  +y«c)2  >T 

°2n+°lQ  <  ’ 

where  Tis  a  fixed  threshold.  Furthermore,  o2nma  +  o2C  very  seldom  becomes  negative, 
allowing  a  modification  to  the  approximate  test 

< 72m  +  crlC  > 

n  a _ a  ^  rj, 

+  *IQ 
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If  all  the  parameters  ma,  <ja,  and  <r„are  known,  this  test  should  yield  near-optimum 
detection.  However,  if  there  is  some  uncertainty,  or  if  the  amplitude  PDF  pa  is  not  really 
Gaussian,  the  major  dependency  of  this  test  is  preserved  by  the  approximation 


where  T is  a  fixed  threshold.  Alternatively,  this  test  can  be  written  as 


M  >  S 

2  rmxm  T,  rm  =  -*=  for  1  <  m  <  M, 

m=l  <  yjQ 


giving  a  local  reference  {rm}  with  no  random  variations. 

In  summary,  for  a  random  signal,  a  near-optimum  processor  is  the  normalized 
crosscorrelator  or  normalized  matched  filter 


M 


M 


r  y  _  m=  1 

m  m  f  M 


—  >  T 

1/2  <  1  3 


regardless  of  the  knowledge  or  statistics  of  scaling  a  (if  positive).  The  only  exception  is  when  a  is 
a  known  constant,  for  which  the  LR  test  is  C  -  jaQ  ^  T\  however,  this  is  not  the  typical  practical 
application. 

The  signal  statistics  are  irrelevant  to  all  these  processor  developments  above,  although 
they  will  affect  the  performances.  This  occurs  because  the  particular  signal  samples,  {sm}  for 
1  <m<M,  are  made  available  to  the  receiver  attempting  to  decide  on  signal  presence  or  absence 
in  that  particular  measurement. 


11/(12  blank) 


PERFORMANCE 


In  general,  the  received  waveform  sequence  is 


=  aasm  +n,  a  = 

m  ffi  m  7 


0  for  Hn 


1  for  Hx  J 


The  average  received  noise  power  is  a2n  -  n2m  (where  the  overbar  represents  an  ensemble 
average),  while  the  received  signal  energy  (when  present)  is  defined  as 

M 

e  =  1L  (asmf=a2Q- 


DETERMINISTIC  TRANSMITTED  SIGNAL 

The  only  processor  of  interest  in  this  subsection  is  the  crosscorrelator  (matched  filter)  with 

output 


MM  M 

C  =  =Z5m(0!aSm  +  nm)  =  «a<2  • 

m= 1  m=  1  m=  1 


For  a  given  value  of  a,  C  is  a  Gaussian  random  variable  with  mean  aaQ  and  variance  Qa2n . 
Therefore,  the  PDF  of  C  is 


pc(u)  =  (2 nQc72)~V2  exp 


{u  -  aaQ)2 

2  Q°2n 


and  the  exceedance  distribution  function  (EDF)  of  C  is 


EM  =  ]dupc(u)  =  ® 


( 


JE 


a - T 

K  on  j 
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where  the  distribution  function  for  the  normalized  Gaussian  variable  is 


<D(x)=  J  dt{ 2n)~V2 


f 


If  a  is  set  at  0  and  then  a  is  set  at  1,  the  false  alarm  and  detection  probabilities  follow  as 


p  =  <D(-  T)  and  Pd=$> 


__  2 

respectively,  for  a  given  value  of  signal  energy  E  =  a  Q. 

If  the  scaling  factor  a  is  a  random  variable,  then  Pd  must  be  averaged  over  a,  using  PDF  pa. 
For  example,  if  a  is  Gaussian  with  mean  ma  and  variance  <r2a ,  then  the  average  Pd  is 


Pd=<D 


4q-t 


i+4 


Q 


(If  oa  -  0,  this  reduces  to  Pd  above,  because  -JE  =  a^Q  -  ma  y[Q  ) 


RANDOM  TRANSMITTED  SIGNAL 

Now,  once  again,  there  are  several  processors  of  interest: 


•  Optimum  processor  C  -  \aQ  <  T 


(known  a). 


•  Normalized  crosscorrelator 


T, 
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•  Crosscorrelator  C  <  J ,  and 

M  > 

•  Energy  detector  Y,xl,  <  T. 

w=  1 


In  order  to  make  progress,  it  is  necessary  to  assign  some  statistics  to  signal  segment  {sm}^ .  It  is 
presumed  that  the  signal  sequence  is  stationary  white  Gaussian,  with  zero  mean  and  variance  a]. 

For  a  given  value  of  attenuation  a,  the  characteristic  functions  (CFs)  of  the  various 
decision  variables  will  be  derived  in  closed  form.  Since  the  average  received  signal  power  is 


of  =  (^m)  =  a1  a) 


the  results  can  be  equally  well  conditioned  on  a2r  or  the  dimensionless  quantity 


R  —  = 


2  2 

a  o. 


which  is  the  received  or  input  signal-to-noise  power  ratio  per  sample.  Thus,  all  the  following 
results  will  be  conditioned  on  a  given  value  of  R.  If  attenuation  a  is  random  and  its  PDF  is 
known,  these  results  can  be  averaged  over  the  statistics  of  attenuation  a  to  find  the  average 
detection  probability. 


The  EDFs  for  all  the  above  processors  cannot  be  found  in  any  convenient  form. 
Therefore,  it  is  necessary  to  resort  to  the  accurate  and  efficient  FFT  procedure  described  in 
the  reference*  for  proceeding  directly  from  the  CF  to  the  EDF.  For  example,  false  alarm 
probabilities  in  the  £-10  range  are  quickly  available  and  accurate. 


*  Albert  H.  Nuttall,  “Accurate  Efficient  Evaluation  of  Cumulative  or  Exceedance  Probability  Distributions  Directly 
from  Characteristic  Functions  ”  NUSC  Technical  Report  7023,  Naval  Underwater  Systems  Center,  New  London,  CT,  1  October 


1983. 
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Optimum  Processor  for  Known  a 

The  decision  variable  for  this  processor  is  given  by 

M  M 

Z  =  C-  \aQ  =  Y,Sm(aaSm  +  nm )  ~  frQ  =  («  “  2 )*Q  +  £  - 

m=l  m=  1 

If  the  signal  sample  values  {sm}  are  held  fixed  for  the  moment,  the  conditional  CF  of  z  is 
/^R,})  =  exp  (i&)  =  exp  [i$(a  - \)aQ  -\^o2nQ]  , 


which  is  based  on  the  fact  that  {nm}  is  a  zero-mean  white  Gaussian  sequence.  Then,  averaging 
over  the  Gaussian  statistics  of  {sm}  results  in  the  CF  of  z  as 

/*(£)'=  I1  +  #2<r>„2  -  i£( 2a  - 1)^2]'^2. 

Since  the  receiver  operating  characteristics  (ROCs)  are  not  affected  by  a  scale  change,  consider 
the  scaled  decision  variable  w  =  z/(crso„),  for  which  the  CF  is 


fM)={  i+f! 


P  =  2a  - 1  = 


- 1  for  Hn 


. 


for  H, 


i  J 


For  a  given  value  ofi?,  EDF  E0(v)  under  H0  is  found  from  f  by  setting  switch  parameter  y9to  -1, 
while  EDF  £i(v)  is  available  by  setting  p  to  1 .  Then,  by  eliminating  threshold  v,  the  ROC,  Pd 
versus  Pf,  can  be  plotted  for  that  particular  value  of  R.  This  procedure  must  be  repeated  for  each 
R  value  of  interest.  (For  R  =  0,  the  two  CFs  are  equal,  leading  to  Pd  =  Pf,  as  expected.) 

An  aid  in  the  numerical  evaluation  of  the  EDFs  discussed  here  is  to  note  that 
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>  Po(“)= a(-w)  , 

leading  to 

oo 

£0(v)  J  #>(*)  =  1  “  ^i(-v)  - 

v 

Thus,  only  Ex(y)  need  be  numerically  evaluated. 

Normalized  Crosscorrelator 

Now,  the  decision  variable  is 


Q  1  M  |  M 

‘WlQ^S-[aaS-  +  n’)=aa^+lQ 


For  any  given  sequence  of  signal  samples  {sm}  for  1  <  m  <  M,  the  last  term  in  z  is  a  Gaussian  RV 
with  zero  mean  and  variance  cr2n ;  thus,  the  conditional  CF  of z  is 


fz{4 |K})  =  exp  (igaa-jQ  -  i£V2) . 


And  since  Q  is  a  sum  of  M  squared  Gaussian  RVs,  its  PDF  is 


uN  1  exp  (-  u!  2a) )  M 

for“>0'  N=~i 


This  leads  to  the  unconditional  CF  of  z  as 


au 

/*(£)  =  exp(-i£2crn2)  J  du  exTp(i^aayfu) 


uN~X  exp  (-«/ 2cr2) 
r(A^)(2c72)W 
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For  scaled  decision  variable  w  =  z/<j„,  and  with  the  change  of  variable  u  =  o]t2 ,  the 
pertinent  CF  is 

fM  =  exP  (-T^2)  |  dt  exp [itajR  t)  -  -f(N)2N-]  ^  ' 

Although  this  result  is  not  in  closed  form,  it  can  easily  and  efficiently  be  numerically  evaluated  for 
any  given  R  by  means  of  a  single  FFT.  Under  H0,  this  reduces  to  fw(%  |//0)  =  exp(-£2/2)  for 
arbitrary  M.  The  alternative  special  case  of  M  =  2  can  be  expressed  as 

im\h i)=exp(-i^)  i+y4f yy  J  • 

Then  Pf  =  <I>(-7),  where  Tis  the  fixed  threshold  against  which  w  is  compared. 

Crosscorrelator 

This  processor  has  decision  variable 

M  M 

z  =  C  =  (aasm  +  nm)  =  ctaQ  +  ^A  • 

m=l  m=l 

The  conditional  CF  of  z  is 

/2(^|R})  =  exp  (i $aaQ  - \?a2nQ) . 

Thus,  using  the  Gaussian  statistics  of  {sm},  the  unconditional  CF  of  z  is 


IB 


With  scaled  decision  variable  w  =  z/(<jsa„),  the  pertinent  CF  is 


As  usual,  a=  0  for  H0  and  a-  1  for#,. 

Energy  Detector 

Finally,  the  last  detector  of  interest  has  decision  variable 

MM  2 

z  =  '£xl='Z(aasm  +nm)  . 

m=\  m=  1 

The  RV  xm  is  zero-mean  Gaussian  with  variance  a a2o2s  +  c2„ .  Using  the  independence  of  the 
{xm}  yields  the  CF  of  z  as 

/I(f)  =  [l-^2(aaV+CT.2)]‘M'1. 

For  scaled  decision  variable  w  =  z/a^ ,  there  follows  CF 

/.(!)  =  [ + 

The  quantity  a  is  equal  to  0  for  H0  and  is  equal  to  1  for  Hh 


For  M  even  and  N  =  M/2,  there  follow  probabilities 


and 


Pd  =  exp 


-T  ) 

<2(1  + P)J 


AT-1  J  /  T  \* 

£  H  V2(l  +  P), 


where  T is  the  threshold  against  which  decision  variable  w  is  compared. 


COMPARISON  OF  DEFLECTIONS 

For  large  numbers  of  samples,  M,  it  can  be  expected  that  the  decision  variables  would  tend 
to  Gaussian.  Then,  if  Pa  =  0.5  were  the  desired  level  of  detectability,  the  deflection 

'  -  ^ 

is  an  important  statistic  in  that  it  tends  to  govern  Pf,  at  least  for  moderate  values  ofP,.  Here,  /4  is 
the  mean  of  the  decision  variable  under  hypothesis  Hk,  while  <r0  is  the  standard  deviation  under  H0. 
However,  since  d  only  involves  first  and  second  moments  of  the  decision  variable,  it  can  yield 
misleading  results  about  which  processors  are  best  in  terms  of  their  ROCs;  the  ROCs  involve  all 
orders  of  moments  and  are  therefore  reliable  measures  for  comparing  processors. 

For  the  optimum  processor,  the  moments  can  be  extracted  from  CF  fw(£)  as 

ff2  =  M1  +  2*)’  d7  =  MYV\R- 

For  the  normalized  crosscorrelator,  N  =  M/2  and 


H  -  a*j2R 


r(N+j) 

T{N) 
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For  large  M,  there  follows  d2  »  MR. 


For  the  crosscorrelator, 

ju  =  aMy[R  ,  a2  =  m{\  +  <z2R)  ,  d2  -  MR . 

Finally,  for  the  energy  detector, 

{i  =  M(l  +  ccR),  a2  =  Ml(\  +  off)2,  d2  =  M\R2 . 

The  deflection  for  the  energy  detector  indicates  a  small-signal  suppression  according  to 
R?',  thus,  for  large  M,  where  small  R  is  of  interest,  the  energy  detector  will  perform  very  poorly 
relative  to  the  other  processors.  This  observation  is  borne  out  by  numerical  evaluation  of  the 
ROCs  of  all  the  processors. 

On  the  other  hand,  for  large  R,  the  energy  detector  has  the  largest  deflection,  which  grows 
as  R2.  Yet,  the  energy  detector  is  not  the  best  processor  for  large  R,  as  the  ROCs  show.  This  is 
an  example  of  the  shortcomings  of  the  deflection  as  a  criterion  for  comparing  processors. 

The  deflection  criteria  for  the  other  three  processors  behave  as  MR  for  small  R,  indicating 
comparable  performance.  However,  for  large  R,  the  deflection  of  the  optimum  processor 
saturates  at  2M,  while  the  deflections  for  the  two  crosscorrelators  increase  as  MR.  Since  the 
optimum  processor  outperforms  both  crosscorrelators,  the  inadequacy  of  the  deflection  is  again 
illustrated. 

In  summary,  the  deflection  criterion  can  be  a  very  misleading  statistic.  It  should  be  used 
only  when  higher  order  moments  are  absolutely  unavailable,  when  the  number  of  samples  is  large, 
when  Pd  «  0.5  is  of  interest,  and  when  Pf  is  moderate  (on  the  order  of  E-2  or  E-A).  In  general, 
comparison  of  the  ROCs  is  the  only  valid  procedure. 
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NORMALIZATION  BEFORE  TRANSMISSION 

When  the  random  signal  segment  {sm}  for  1  <  m  <  A/ is  drawn,  suppose  that  it  is 
normalized  prior  to  transmission,  according  to 


for  1  <m<  M  , 


where  ET is  the  transmitted  signal  energy.  Since  the  receiver  knows  {sm},  it  has  no  uncertainty 
about  the  transmitted  signal,  and  this  case  reverts  to  the  deterministic  signal,  even  though  the 
particular  transmitted  sequence  varies  from  trial  to  trial.  That  is,  the  transmitted  deterministic 
signal  varies  from  trial  to  trial. 

The  received  signal  waveform  under  Hx  is  asm,  where  a  is  the  medium  attenuation.  The 
LR  test  is,  as  earlier. 


c<r  , 


and  the  performance  is  again 


P.  =  fl>(-  T) ,  Pd  =  <t 


J 


for  a  given  value  of  the  received  signal  energy 


M  2 

£  =  Z(aSm)  =a2Er. 

m=l 
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GRAPHICAL  RESULTS 


The  ROCs  for  detection  of  a  deterministic  signal  are  given  in  figure  1.  Because  the 
received  signal  is  not  random  in  this  case,  these  curves  are  parameterized  in  terms  of  the  total 
received  signal  energy  E.  In  particular,  a  decibel  (dB)  measure  relative  to  the  average  received 
noise  power  is  utilized.  These  characteristics  are  parallel  straight  lines  for  all  decibel  values. 

All  the  remaining  ROCs  are  for  random  transmitted  signals  and  are  parameterized  with 
respect  to  the  received  average  power  signal-to-noise  ratio  R,  namely,  i?(dB)  =10  logi0(/?).  The 
situation  of  a  low  number  of  samples,  M=  10,  is  depicted  in  figure  2  for  the  optimum  processor. 
(The  quantity  Aj,  is  the  increment  in  £  used  in  computing  the  EDF  directly  from  the  corresponding 
CF;  Mb  and  b  are  shift  values  guaranteeing  a  positive  random  variable;  and  Mf  is  the  size  of  the 
FFT  employed.  These  computational  parameter  values  are  presented  in  case  the  user  desires  to 
confirm  or  extend  these  results.)  The  corresponding  results  for  the  optimum  processor  with 
M=  100  and  M=  1000  are  displayed  in  figures  3  and  4,  respectively,  along  with  the  required 
values  of  the  computational  parameters  that  guarantee  accurate  false  alarm  and  detection 
probabilities. 

Similar  ROCs  are  presented  in  figures  5,  6,  and  7  for  the  normalized  crosscorrelator  with 
M =  10,  100,  and  1000,  respectively.  The  corresponding  results  for  the  (unnormalized)  cross¬ 
correlator  are  given  in  figures  8,  9,  and  10,  while  the  energy  detector  results  are  depicted  in 
figures  11,  12,  and  13. 

A  low-quality  operating  point  is  defined  as  Pf  =  \E-3  and  Pd  =  0.5,  while  a  high-quality 
operating  point  is  defined  as  Pf  =  IE-6  and  Pd  =  0.9.  The  first  observation  to  make  is  that  the  two 
“usual”  processors,  the  energy  detector  (ED)  and  the  crosscorrelator  (CC),  cannot  be 
unequivocally  ordered  in  terms  of  performance.  For  example,  for  the  small  time-bandwidth 
product  of  M=  10,  reference  to  figures  8  and  1 1  shows  that  the  CC  outperforms  the  ED  at  the 
low-quality  point,  requiring  a  signal-to-noise  ratio  per  sample  of  1 .6  dB  versus  3.3  dB.  On  the 
other  hand,  at  the  high-quality  point  for  M=  10,  the  ED  is  better  than  the  CC,  requiring  only 
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9.3  dB  versus  12.3  dB.  However,  for  larger  numbers  of  samples  M,  such  as  M-  100  in  figures  9 
and  12,  orM-  1000  in  figures  10  and  13,  the  CC  considerably  outperforms  the  ED  for  all  signal- 
to-noise  ratios. 

Tables  1  and  2  below  list  the  values  of  the  input  signal-to-noise  ratio  per  sample  (in 
decibels)  required  in  order  to  achieve  the  low-quality  and  high-quality  performance  levels  for  the 
ED  and  CC,  as  well  as  for  the  normalized  crosscorrelator  (NC)  and  the  optimum  processor  (OP). 
It  can  be  seen  from  these  results  (or  figures  2,  3,  and  4  relative  to  figures  5,  6,  and  7,  respec¬ 
tively)  that  the  NC  performs  virtually  at  the  optimum  performance  level  attainable  by  the  OP.  The 
only  discrepancy  is  0.2  dB  for  M-  1 0  at  the  high-quality  operating  point.  This  is  particularly 
striking  when  it  is  recalled  that  the  NC  does  not  need  any  additional  information  about  the 
medium  attenuation,  while  the  OP  does  need  and  use  this  information. 

The  rightmost  column  in  the  tables  is  derived  from  the  deterministic  signal  results  in 
figure  1.  It  does  not  constitute  a  strictly  legal  comparison  case;  it  has  been  added  only  for 
completeness  and  curiosity.  Furthermore,  the  values  tabulated  have  been  converted  to  a  per- 
sample  value  by  subtracting  10  logio (M)  from  the  values  obtained  in  figure  1,  which  measure  the 
total  received  signal  energy  E.  Nevertheless,  the  values  in  this  last  column  are  only  slightly 
smaller  than  those  for  the  OP  and  the  NC,  indicating  rather  close  performance  levels  for  the 
deterministic  versus  random  signal  cases. 


Table  1.  Required  Signal-to-Noise  Ratio  for  P/=  IE-3  and  Pj  =  0.5 


M 

ED 

CC 

NC 

OP 

Deterministic 

10 

3.3 

1.6 

0.0 

0.0 

-0.2 

100 

-3.0 

-10.0 

-10.2 

-10.2 

-10.2 

1000 

-8.4 

-20.2 

-20.2 

-20.2 

-20.2 
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Table  2.  Required  Signal-to-Noise  Ratio  for  Pf  =  IE-6  and  Pi  =  0.9 


M 

ED 

CC 

NC 

OP 

Deterministic 

10 

9.3 

12.3 

7.3 

7.1 

5.6 

100 

0.8 

-3.5 

-4.2 

-4.2 

-4.4 

1000 

-5.2 

-14.3 

-14.3 

-14.3 

-14.4 
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Probability  of  False  Alarm 


Figure  1.  ROCs  for  Deterministic  Signal 


Probability  of  Detection 


30 


0.999 


Probability  of  Detection 
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512 


0.0025 


512 


SUMMARY 


The  characteristic  functions  of  the  decision  variables  for  the  four  different  processors  used 
to  detect  random  signal  transmissions  have  been  derived  in  closed  form  and  then  numerically 
exercised  to  yield  the  corresponding  accurate  receiver  operating  characteristics.  The  number  of 
data  samples  M  available  during  the  observation  interval  is  arbitraiy,  as  is  the  input  signal-to-noise 
ratio. 


The  usual  energy  detector  and  crosscorrelator  have  operating  characteristics  that  cross 
each  other,  the  exact  locations  depending  on  the  value  of  M  and  the  particular  signal-to-noise 
ratio.  Neither  of  these  processors  performs  as  well  as  the  normalized  crosscorrelator,  which  itself 
performs  within  a  small  fraction  of  a  decibel  of  the  optimum  processor.  The  obvious  conclusion  is 
to  recommend  use  of  the  normalized  crosscorrelator  as  the  preferred  detector,  at  least  under  the 
conditions  utilized  in  this  study. 
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